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We present an analytical model describing the full electromagnetic propagation in a THz time-domain 
spectroscopy (THz-TDS) system, from the THz pulses via Optical Rectification to the detection via Electro 
Optic-Sampling. While several investigations deal singularly with the many elements that constitute a 
THz-TDS, in our work we pay particular attention to the modelling of the time-frequency behaviour of all 
the stages which compose the experimental set-up. Therefore, our model considers the following main 
aspects: (i) pump beam focusing into the generation crystal; (ii) phase-matching inside both the generation 
and detection crystals; (iii) chromatic dispersion and absorption inside the crystals; (iv) Fabry-Perot effect; 
(v) diffraction outside, i.e. along the propagation, (vi) focalization and overlapping between THz and probe 
beams, (vii) electro-optic sampling. In order to validate our model, we report on the comparison between the 
simulations and the experimental data obtained from the same set-up, showing their good agreement. 

Generation and detection of THz electromagnetic waves (0.1-10 X 10 12 Hz) is not a novel topic 1 " 3 , but it 
has been arousing an ever-increasing interest only in the last decade. Recent studies demonstrated that 
many substances and particularly bio-polymers like proteins, amino and DNA possess specific global and 
sub-global modes in the THz band. Hence, particular consideration has been recently devoted to THz spectro- 
scopy since it reveals information on the conformational stage of molecules and potentially enables their dis- 
crimination in various compounds. Indeed, THz pulses have been broadly used, ranging from spectroscopy and 
time-resolved pump-probe experiments to imaging applications 4 ' 5 . Moreover, the negligible ionization power of 
the THz radiations compared to optical and X-rays technologies make them perfectly suitable for biological 
applications 6 . For the above-mentioned reasons, THz technology are penetrating in lots of areas, beyond fun- 
damental Physics investigation, spanning from medical care and homeland security to cultural heritage conser- 
vation 7 " 13 . The most challenging technology development remains in the area of THz sources and detectors, since, 
in general, both wide band and high power are required at the same time. Thanks to the improvement of the 
pulsed laser technology, different types of stable THz sources have been developed so far, with peak power values 
ranging from kW to MW 1415 . In particular, Optical Rectification (OR) 16 combined with Electro-Optic Sampling 
(EOS) detection is still the configuration of choice in the few standard high energy pulsed THz system. For this 
reason, in this paper we present an analytical model able to simulate the full electromagnetic propagation which 
takes place in the THz-TDS set-up 17 , sketched in Fig. 1 (refer to Methods), based on OR and EOS and employing 
<110> ZnTe crystals. We will explain and justify the approach and the approximations which underlie the 
modelling of each stage and that allowed us to build a complete transfer function in the frequency domain, in 
order to predict the set-up behaviour. 

Results 

Optical rectification. In noncentrosymmetric crystals the second- order nonlinear optical susceptibility x (2) is 
nonzero 18 and, for this reason, parametric processes (such as second-harmonic generation, sum- and difference- 
frequency generation, parametric solitons, optical rectification, etc.) can take place 19 " 24 . As regards Optical 
Rectification, when a CW laser impinges on such a material, the latter responds generating a DC -i.e. 
rectified- electric field proportional to the pump intensity. On the other hand, when an ultrashort (i.e., 
broadband) laser pulse (—100 fs) passes through a non-centrosymmetric crystal, optical rectification performs 
a mixing of the frequencies within each pulse -i.e., difference-frequency generation (DFG) between pairs of 
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Figure 1 | Sketch of a typical THz Time Domain Spectroscopy (THz- 
TDS) set-up. 

frequencies within the pulse bandwidth- generating a low-frequency 
radiation, that is the THz wave. The time dependent equation 
describing this phenomenon is 25 " 27 : 



c 2 dt 2 



e 0 c 2 dt 2 



where E THz is the instantaneous THz field, P THz is the linear 
polarization induced by the THz electric field and P^ z is the 
nonlinear polarization induced by the electric field of the laser 

pump Epp(f,t) = 1/ 2£(r,t) exp ^ {co PP t — k PP -r^j j + ex., where w PP 
and k PP are the carrier frequency and the wave vector respectively, 
and £(r,f) is the slowly- varying envelope both in r and t, for each 
pulse. We will neglect the contribution of other nonlinear processes, 
such as, second harmonic generation, two photon absorption, free 
carrier absorption and self-phase modulation. In addition, the pump 
depletion is always negligible in most practical cases, since the 
conversion efficiency for the OR case (~10~ 6 -10~ 4 ) is far below 
the upper limit of the quantum efficiency, given by the ratio of the 
central THz frequency and the optical carrier. By assuming a near- 
infrared pumping and pulses of transform-limited duration of 100 fs 
or longer, we neglect the second time -derivative of £ (slowly varying 
envelope approximation, SVEA). Furthermore, we consider both 
pump and THz beams as plane waves, neglecting the dependence 
on the transverse coordinates (x, y), being the diffraction length in 
those cases always much larger than the crystal thickness in all the 
generated spectrum. By taking the Fourier transform of eq. (1), with 
P L THz (Q) = 8oX (1 \^)Ethz(zM we obtain 28 : 



dz 2 + c 2 



E THz (z,Q) -- 



+ 00 

v(2)/YVk f & 



d z \E PP (t) \ 
dt 2 



2 

low n 
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where the over-lines stand for the Fourier operator and 

nmz(Q) = a/1 + is me refractive index at the THz frequency 

(Q). It is worth highlighting the presence of the second spatial derivative 
in eq. (2), since the generated THz wave has a bandwidth comparable to 
its central frequency, i.e. SVEA does not hold for THz radiation. 
According to the laser features commonly employed in the THz 
collinear generation, we approximate the pump pulse envelope with a 

Gaussian profile \E PP {z,t)f hw = £q/2 exp — (t — z/v g ) 2 / T pp] > where 

by definition t pp = t P whm /2 y/ln (2), t fwhm being the full-width at 
half maximum pulse duration of the pump intensity, v g = cln g 
and Eq the group velocity and the magnitude of the electric field, 
respectively. Here, we assume that t pp is not varying along the 
propagation because of dispersion, since the thickness of the 



crystal (~1 mm) is much smaller than the dispersion length of 
both fields. Thus, by considering the boundary conditions 
E THz (0,Q) = dETHz(.z,Q)/dz\ z=0 = 0, one obtains the exact solution 
for the THz electric field in the frequency domain, from eq. (2): 



E THz (z,Q) 



Elx^Tpp^/n 
(n 2 THz {Q)-nl) 



exp 



(3) 



n T Hz(ty 



+ 2 u + 



n T Hz(ty 



THz wave properties. In eq. (3), we wrote the THz frequency- 
dependent electric field as the product of three main terms in 
order to point out the properties of the THz wave. The third factor 
shows that the THz waveform is actually the result of the interference 
among three waves: the first two travel with the same phase velocity 
VTHz(Q) = c/Re{nTHz(Q)} (precisely, one is backward 29 ' 30 and the 
other one is forward with respect to the z axis), whereas the third 
one is a forward wave with phase velocity matching the group 
velocity of the pump beam. However, the amplitudes of the three 
contributions are quite different: it results that although the 
backward wave vanishes at the frequency satisfying the phase- 
matching condition, as explained below, it has always negligible 
magnitude in the whole THz range. The relative term will be 
therefore neglected. Conversely, the forward wave reaches its 
maximum at the phase-matching. It is generally comparable with 
the third wave in the whole THz range and they together contribute 
to the generated THz beam. The second factor in eq. (3) is a Gaussian 
term centred at Q = 0, which acts as a low-pass filter: the shorter the 
transform-limited pulse duration T pp the wider the bandwidth. 
Indeed, there is an upper limit for the THz band Av THz obtainable 
via ultrashort laser pulses, according to Heisenberg's principle 25 , i.e., 
the time-bandwidth product, Avthz = (2 In 2)/ (tctpp)— 0.441 /t pp . 
However, in most practical cases of femtosecond-pulse driven 
generation, the achievable bandwidth is never as wide as it is 
determined by the latter expression, since it mostly depends on the 
chosen crystal, as we will show later. Finally, the first factor in eq. (3) 
is comprehensive of the nonlinear parameter and the difference 
between the refractive index of the crystal in the THz range and 
the infrared group index. Both x (2) (^) and wthz(^) = n(Q) + w(Q) 
are complex frequency- dependent functions, accounting for dis- 
persion and absorption, whose effects are not negligible because of 
the ultrashort nature of the pulse pump. Hence, we model the 
refractive index of the crystal, in the THz range, using Maxwell- 
Helmholtz-Drude's semi-quantic model 31 of the transversal phon- 
on resonance modes which establishes that the dielectric permittivity 
can be expressed by: 



WO) = £ (oo) + If^% 
il T0 — il +iy TO il 



(4) 



where s(oo) is the dielectric constant at frequencies much higher than 
1 THz, e(0) is the DC dielectric constant, Q TO is the transversal 
optical (TO) resonance frequency, y TO is the full-width of the 
mode. Noteworthy, this model is the best choice for our purposes, 
since ZnTe is featured by a very weak coupling between the 
transversal and longitudinal modes. Hence, for frequencies close to 
Q ro > me resonant absorption becomes so strong that the material 
appears completely opaque. As Q TO /2n = 5.32 THz 32 , ZnTe sources 
exhibit a contiguous generation window between hundreds of 
GHz up to 3 THz. Note that other materials guarantee larger 
bandwidths 33 " 36 . ZnTe is still the material of choice in many set-up 
due to a nicely shaped THz generation spectrum and the significant 
power-spectral density within the lower THz band. For the nonlinear 
susceptibility, we have used Miller's rule 18 . 
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Figure 2 | (a) Coherence Length vs. frequency for the case of three different pump wavelengths in ZnTe. The black line indicates the generation crystal 
thickness, (b) Sketch of the emitting cone for the THz beam both inside and outside the ZnTe crystal. 



Phase-matching and coherence length. The phase matching 
condition is fulfilled when the optical group velocity matches the 
THz phase velocity: 



vthz(&*) = vJX*)on(Q* = 2nv*) = n g (X*) 



(5) 



For ZnTe it turns out that n*„ = n% 



-3.285, v* « 2.55 THz, 2* 



788 nm 37 . The THz frequency O* is generated with the maximum 
possible conversion efficiency. This could be further ascertained by 
manipulating eq. (3), as it follows: 



Eth z (z,Q)ccL, 



'gen 



sin(Ak(Q)L gen /2) 



Ak(Q)L gen /2 



(6) 



is the 



where Afc(O) = Q(n — n g )lc is the mismatch term and L gen 
generation crystal thickness. Close to the PM condition, i.e. A/c(0) ~ 
0, the THz field magnitude scales almost proportionally to the crystal 
length. Anyway, for each frequency within the bandwidth and for a 
fixed laser wavelength, there is an upper theoretical limit for L gen , 
related to the so-called coherence length: 

A THz 



L co h(Q) = 



2\n g -n(Q)\ 



(?) 



The physical meaning of (7) is that, when the propagation length 
approaches L coh , the generation process equates the back- conversion 
process that reconverts the THz photons into pump power until the 
THz field vanishes for z = 2L coh . This means that it is necessary to 
achieve a trade-off between power and bandwidth when the 
thickness of the generation crystal has to be chosen according to 
the proper application. To this end, Fig. 2a shows how the 
coherence length is related to the frequency for three different 
pump wavelengths. We note that employing a 1-mm thick ZnTe 
crystal, the bandwidth is maximum (about 2.9 THz) for a 778 nm 
pump (red line), but it shrinks when the laser wavelength moves 
away from A*, even if the generation efficiency (related to L coh ) 
becomes higher for all the frequencies within the band (blue and 
green lines). 

The L gen = 2 mm case is also plotted showing that the generation 
process is not efficient for a relatively large THz window when 2* is 
employed as pump source. Indeed, in this case, as pointed out by the 
double brown dashed-arrow, the most of the bandwidth achievable 
with ZnTe in the rage O < O ro > is not practically generated, resulting 
in a narrow-band spectrum centred at the frequency O*. 

Crystal interface. The effect of the crystal-air interface can be 
formalized considering the dispersion of the Fresnel reflection and 



transmission coefficients within the THz band within the hypothesis 
of non-diffracting (collimated) beam. These considerations are 
sketched on Fig. 2b. We can then introduce frequency- dependent 
Fresnel coefficients for transmission (T) and reflection (R), as: 



T^ 2 (ft): 



2y/firffe(fl) 

y/8TH Z (Q) + l 9 



y/firffe(fl)-l 

y/8TH Z (0) + l 



(8) 



where "1" and "2" in the subscripts stand for ZnTe and air, 
respectively. As the dielectric constant diverges for O — > O ro , eq. 
(8) has a resonant complex nature that emphasizes reflection in the 
high THz band. From the previous consideration, it immediately 
follows that an important phenomenon which takes place inside 
the crystal is the Fabry-Perot (FP) effect, which acts as a selective 
filter, according to: 



Tpp(Lgen,Q) 



T^ 2 (Q)exp[ 


-ai»nfe(0)5= 




l-i^ 2 (Q)exp 




gen 
C J 



(9) 



From eq. (9) we note that if the crystal is thin enough, the 
longitudinal resonances are established for those frequencies which 
do not fall into the achievable THz bandwidth, thus resulting in no 
THz pulse-repetitions outside the material. In addition, eq. (9) 
somewhat modifies the THz bandwidth, because it mostly 
transmits those frequencies which fulfil the nonlinear condition, 
2Qn(Q)L gen /c + ^(O) = 2 nm, where 0(0) is the phase of 
£ 2 (0)and m is the order of the longitudinal mode. 

Diffraction. Outside the crystal, the broadband THz beam pro- 
pagates toward the mirror that lies in the far field region, but it 
diffracts out much quickly than a Gaussian beam. Therefore, we 
cannot use Rayleigh's law to quantify its aperture along the pro- 
pagation, since the latter only holds in the paraxial approximation. 
Instead, we can apply Bethe's theory of diffraction by small holes 38 . If 
fact, since the THz field is proportional to the pump intensity, we can 
assume that the THz wave is generated by a spot size whose waist 
satisfies the following relation: Wthz = w pp I V2. This means that the 
THz source size could be smaller than most of the emitted wave- 
lengths, as Bethe's theory assumes. Avoiding the exact treatment of 
the diffraction, the main result of the theory is that we can introduce a 
different confocal parameter similar to the Rayleigh length, defined 
as z B (Q) = k THz w TH j2 = Q 2 w PP /(4V2c 2 ), which shows how the 
lower frequencies diffract quicker than what is predicted by 
Rayleigh's law, since in this case z B depends on the square of 
frequency. However, we note how for larger waist size w PP , Bethe's 
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Figure 3 | (a) I T Di ff\ 2 transfer function (eq. 9) compared with Rayleigh's and Bethe's diffraction (plotted for a pump beam waist fixed at w PP ~ 30 fim). 
Sketch of the (b) collimating stage and (c) focusing stage. 



diffraction is practically negligible for the whole bandwidth, as we 
will show in the following. Then, in order to compute the field rate 
gathered by the first mirror, we would follow the procedure proposed 
in Ref. 28, calculating the transfer function for the diffraction stage: 



T Di ff(Q,6 mir ) = 



\ 



1 — exp 



tan 2 (0„ 



W THz 



(10) 



here 0 mir = arctan[l/(2/#)] and/ # =f/D are the semi-aperture and the 
/-number of the mirror, respectively, whose diameter and focal 
length are D and /, respectively (see Figs. 3b and 3c). For more 
accuracy, we distinguished two cases, defining the confocal 
parameter z di ^Sl) = z B if k THz w THz < 1 whereas z^Q) = z R if 
krHz w THz > 1> with Zr = 1cthzWthz/2 me usual Rayleigh length, 
which anyway describes quite well the behaviour of the higher 
frequencies. In order to show the effect of Bethe's diffraction, we 
report on a plot of eq. (10) in Fig. 3a, under our working 
conditions, i.e. f# ~ 1.4, w pp ~ 30 fim, which shows the different 
behaviour of the two laws. As we can see, the frequencies lower than 
2 THz -which fall in the so-called sub -wavelength regime- are very 
weakly collected by the mirror since its numerical aperture is not 
wide enough. The black circles in Fig. 3a indicates the trend switch of 
the curve, i.e. those frequency whose propagation follow again 
Rayleigh's law. Since we are mainly interested in as diffraction is 
related to the pump beam waist and how the former affects the 
THz pulse bandwidth, we simulated several cases, as reported in 
Fig. 4, where the generation parameters are fixed, except for the 
pump waist size. In particular, in Fig. 4a, the comparison among 
the bandwidths collected through the first parabolic mirror - 
varying the parameter w pp - is shown. Again, it is possible to 
recognize how the mirror cannot efficiently collect the lowest side 
of the THz bandwidth - which includes the sub -wavelength regime, 
highlighted by the dashed lines- and that smaller w PP , wider such 
region, as it is stressed by the trend of the curves. The black thicker 
curve indicates our generation condition, meaning that most of the 
bandwidth we generated falls into Bethe's regime, so its inclusion in 
the model is a crucial aspect in order to achieve agreement with the 
experiments. Simulations reasonably show that the generation 
process -for a fixed pump power- is more efficient when the pump 
beam is tightly focused (in fact, the THz intensity for w PP > 200 fim 
is at least one order- of-magnitude lower than the other cases). 
However, it should be kept in mind that waist sizes too small (w PP 
< 200 fim) fail to satisfy the hypothesis we imposed when we solved 
eq. (2), and which could probably give rise to other nonlinear 
phenomena due to the extreme pump beam focusing, phenomena 



which are not included in the model and which we are not interested 
to. On the other hand, pump waists too large (w PP > 500 fim) are not 
suitable -in addition to the too low efficiency conversion- because, for 
this case, it has been shown that the THz radiation generated at 
peripheral transversal positions along the pump beam destructively 
interfere with the one generated on the axis 39 . It is also worth noting 
that the THz pulse gathered by the off- axis is featured by a central 
frequency (given by the intensity- weighted average frequency, i.e. the 
first- order moment with respect to intensity) that moves towards 
higher frequencies as the waist size decreases. 

Focusing. THz beam focusing is hampered by the same issues of the 
collecting stage. Because of the broadband nature of the THz beam, 
each frequency- component is being focused not only to different 
spot sizes, but also to different focal planes. This is a serious prob- 
lem because the strength of the THz electric field is distributed on a 
volume of the detection crystal, making the detection process more 
difficult, being based on the optical index modulation caused by the 
THz field, as we will see later. However, we have to stress that the 
bandwidth of the beam arriving at the second mirror has not 
significant low-frequency terms, so their contribution is marginal. 
This consideration allows us to simply apply the geometrical optics 
rules. The beam size arriving at the second mirror is: 



»(«)= 



D 



2^2' 



p 



(11) 



where we have supposed, for sake of generality, that the THz beam 
could diffract more than the dimension of the mirror. Actually, eq. 
(11), plotted in Fig. 4b, shows how in the sub-wavelength regime, the 
THz beam diffraction is so strong that the beam radius - ideally 
measured in front of the mirror - is larger than the cross section of 
the mirror. This means that the collected THz beam waist is 
essentially fixed by the mirror dimension for most of the band of 
interest, for several suitable values of w PP . Now, if we use Gauss' 
formula, it is found out that the mirror focuses the THz beam to a 
waist as wide as ry 0C (Q) = 2cf/Qr in (Q), from which the focusing effect 
could be modelled as 40 : 



rjoM _ 2cf 



r in (Q) Qff n (Q) 



(12) 



A plot of this transfer function is reported in Fig. 4c, revealing some 
interesting aspects. Certainly, the most important effect of this stage 
is the low- frequency enhancement of the generated bandwidth, as it 
can be also seen in Fig. 4d, in the case of r in practically constant within 
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Figure 4 | (a) Comparison among the bandwidths collected through the first off- axis for different values of the pump beam waist w pp (normalized with 
respect to the maximum of the curve at w pp = 10 (am). It is also shown Bethe's regime extension (dashed arrows), (b) Frequency-dependent transversal 
size of the THz beam collected through the first off-axis (normalized with respect to the effective mirror radius R e ff =D e ff/2 = D/(2y/2)). Curves 
corresponding to a waist smaller than 50 (am are constantly equal to 1 for the whole investigated THz range, (c) Transfer function I Tf oc \ 2 of the focusing 
stage through the second off-axis. Curves corresponding to a waist smaller than 50 um are entirely overlapped, (d) Comparison among the bandwidths 
focused through the second off- axis for different values of the pump beam waist w pp . Each curve is normalized with respect to its own maximum. In all the 
figures, the thicker black curves highlight our working conditions. 



the entire band. This because the focusing effect is stronger for wide 
rather than thin beams impinging on the mirror. However, we want 
to highlight that, for larger pump waist (w PP > 200 fim), r in is no 
longer constant within the whole band but tends to decrease starting 
from that frequency whose beam component does not diffract out 
more than the mirror size. Considering eq. (12), this entails that the 
higher frequencies are now emphasized as well (see, for example, the 
Wpp = 500 fim case in Fig. 4d). Consequently, this stage actually 
performs a reshape of the spectrum such that, for example, the 
normalized bandwidth for the w PP = 500 fim case is almost over- 
lapped with those for the w PP < 30 \im case (even if the spectra 
amplitudes are obviously more than one order- of- magnitude 
different). Anyway, we note how the most interesting aspect is that 
if it is possible to reach a trade-off for the pump waist size (as the w PP 
= 300 /urn case), the reshape-effect of this stage could be helpful to 
achieve wider bandwidths for fixed nonlinear generation crystal. 

Electro-optic sampling. Free-space electro-optic sampling is a 
coherent pulse detection technique, which has been extensively 
experimented for several years 41 " 43 . The working principle is sket- 
ched in Fig. 5a (refer to Methods). The advantage of EOS is that one 
can determine not only the amplitude of the pulses but also their 
phases with high accuracy (< 10" 2 rad). Although the simplicity of 
the above-described mechanism, EOS is afflicted by the same issues 
of the OR, such as 44,45 : (i) dispersion and absorption inside the crystal 
and at the air-crystal interface; (ii) Fabry-Perot effect inside the 
detection crystal; (iii) finite pulse duration of the probe beam; (iv) 
phase-mismatch between probe and THz pulses. Furthermore, a 



critical problem is (v) the geometrical overlapping between the 
probe and THz beams. While the modelling of (i-ii) is the same as 
the previous stages, for (iii-iv) we write the following frequency 
response for the detection stage: 

F(Q,co PB ) =A opt (Q)x (2 \co PB] n,CD PB - ft)A0> (Q) (13) 

where A opt (Q)is the autocorrelation of the Gaussian probe pulse, 
which accounts for the loss of resolution due to the finite pulse 
duration, x {2) ( w pb> Q w pb — Q) is the second-order susceptibility 
for the DFG case (the same as the OR case, according to Kleinman's 
symmetry) and A<I>(Q) is the frequency- dependent mismatch factor. 
With simple manipulations, it is found out: 



\Ztz 

A opt (Q) = — exp 



Q 2 T 2 



AO(Q) : 



exp ( - iAk(Q)L det ) - 1 
iAk(Q) 



(14) 



(15) 



where t pb is the probe pulse duration, A/c(Q) = k P (w PB ) + k(Q) — 
kp(co PB — Q) is the mismatch term between the probe field, the THz 
field and the frequency- difference field wave vectors, respectively, 
and L det is the thickness of the detection crystal. 

Beams overlapping. The issue mentioned in (v) is due to the fact that 
THz focusing is frequency- dependent as well. The THz beam could 
be ideally considered as the overlap among lots of single-frequency 
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Figure 5 | (a) EOS working principle for a train of pulses, (b) Sketch of the overlapping between probe and THz beams inside the detection crystal. 



beams, whose waist is frequency- dependent. So, those frequencies 
whose spot size is more spatially overlapped with the probe one are 
being detected with higher efficiency, as shown in Fig. 5b. To 
compute this effect, we have to consider the transversal beam 
profiles and solve the overlap integral: 



T over i ap (Q)cc \E PB (r)\ 2 E THz (r)dr = 



\w 2 PB + 2rj 0C (Q) 



(16) 



where we have assumed a Gaussian shape for both beams, whereas 
rf 0C (Q)is the THz waist and w PB is the probe one. It should be noted 
that eq. (16) is only valid when the two beams are exactly concentric. 
This effect can be quite important, since the probe beam size is 
generally smaller than the THz one, due to its focusing into the 
detection crystal, in order to obtain an enhancement of the EOS 
efficiency. 



Discussion 

As we mentioned in the beginning, we can then multiply together all 
the previous terms in order to obtain an overall transfer function 
T 

1 set-up 



in the frequency domain, which models the entire set-up 



set — up(Q>f#>@mir)L(i e t ,(Dp B ,Tp B ) 
X E g THz(QL g en,ti>pp,Tpp) 



(17) 



behaviour providing the detected THz electric field: 



where Ej^ z is the THz electric field emerging from the generation 
crystal. Obviously, taking the Fourier transform of eq. (17), one can 
analyse the instantaneous THz electric field in the temporal domain. 
However, it should be kept in mind that this can be only performed 
with a numerical approach, since generally E^ z (t) cannot be 
expressed with elementary functions. In order to test the effective- 
ness of the model, we report on the comparison with the experi- 
mental results. As we have mentioned through the text, we built an 
ad-hoc model for our set-up, so we have already given many experi- 
mental parameters, which we list again in the following for simpli- 
city. The generation and detection crystals are both <110> ZnTe: 
the first one length is L gen = 1 mm, whereas for the second one L det = 
2 mm. The optical source is a Ti:Sapphire pulsed laser (Mai Tai, 
Spectra- Physics) whose centre emission line is k PP = 800 nm. The 
intensity pulse duration is t pp = t pb = 200 fs, whereas the waist sizes 
focused to the crystals are equal to w PP = w PB = 30 fim, for both 
beams. The/- number of both collimating and focusing optics isf# ~ 
1.4. By using these parameters for the simulations, we obtained the 
results plotted in Figs. 6a and 6b, which represent the comparison 
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between the experimental (circles) and the simulated (red solid-line) 
spectra and waveforms, respectively. The experimental data were 
obtained in a temporal window of about 50 ps and are indicated as 
green circles in both the figures. In order to make a convenient 
comparison with the model, we had to perform a post-processing 
of these data. In fact, we note how the corresponding spectrum in 
Fig. 6a (green circles) is corrupted by some unwanted phenomena, 
mainly related to the choice of the observation window and to the 
laboratory atmosphere. As regarding to the first point, the obser- 
vation window at the detection stage has to be wide enough in order 
to achieve a good frequency resolution. However, in this way, the 
instrumentation also reveals the spurious THz pulse repetitions, 
because of the FP effect inside the crystals. From theory, we know 
that the Fourier transform of a periodically repeated signal in time 
domain, appears as the sampled spectrum of the fundamental wave- 
form (which, instead, we are interested to). With regards to the 
standard laboratory atmosphere, it results quite opaque to THz 41 , 
because of the presence of many resonance lines in such a range of 
the spectrum, due to the water vapour. As a consequence, the band- 
width lacks the contribution of such frequency components, because 
of the efficient resonant absorptions. In particular, we stress that the 
notch around 1.7 THz in Fig. 6a is so strong that the waveform is 
featured by a residual quasi-sinusoidal ripple at that frequency for 
positive times. Obviously, the latter is not predictable with an ana- 
lytical approach. These considerations help us to explain the form of 
the curve (indicated by the green circles). Nevertheless, to test our 
model, we needed to manipulate a quasi-ideal waveform and its 
related spectrum. Thus, we isolated the fundamental THz pulse in 
time domain - corresponding to the blue circles in Fig. 6b - and 
calculated the corresponding spectrum plotted via the blue circles, 
as shown in Fig. 6a. The two curves can now be compared with our 
simulations. As we can see, the agreement between the model (red 
line) and the experimental curves is very good, meaning that our 
model is able to describe quite well the various phenomena that take 
place in the set-up. In particular, the spectra demonstrate the hypo- 
thesis of sub -wavelength generation for almost all the bandwidth, 
exactly as described by Bethe's diffraction, which is a crucial aspect of 
the THz propagation modelling. Moreover, it can be noticed that the 
central frequency is about f mean = 1 THz, whereas the bandwidth 
measures about B 1/e = 1.2 THz, as provided by the PM conditions 
(5) and (7) with L gen = 1 mm. Conversely, in time domain, the 
matching between the simulated and experimental waveforms is very 
good, and both show that the obtained THz profile significantly 
deviates from the single-cycle case -as usual for ZnTe crystal- since 
the bandwidth generated is not wide enough. 

In this work we built an analytical 1-D model able to simulate a 
typical generation and detection THz pulses set-up based on OR and 
EOS. We studied the phenomena that occur along the THz beam 
propagation, in order to obtain an analytical approach that allowed 
us to predict the behaviour of each stage by means of a proper 
frequency- dependent transfer function. We compared the results 
of the simulations with the experimental data, obtaining a good 
agreement. This means that our study of the theoretical phenomena 
was careful enough to take into account the main frequency prop- 
erties of the stages. The peculiarity of the simulated waveforms con- 
sists in assuming them as reference signal, when employed for the 
study of particular properties of materials such as their frequency- 
response (for example in the TDS technique), whenever it is not 
possible to obtain their profile without external and unwanted phe- 
nomena. 

Methods 

Experimental set-up. The experimental set-up is sketched in Fig. 1. A train of 
ultrashort laser pulses at 800 nm, emitted by a Ti:Sapphire laser (Mai Tai, Spectra- 
Physics) is focused through a lens into a <110> ZnTe sample, which acts as the 
nonlinear generation crystal, giving rise to the OR. It is worth to highlight the 
possibility of tailoring the size of the pump waist focused into the generation crystal, 



by moving the position of the lens along the propagation axis. The resulting THz 
broadband beam and the optical beam emerge from the crystal and are handled by a f- 
2f-f coherent system consisting of two 90° off- axis parabolic metal mirrors with focal f 
= 50 mm and diameter D = 50.8 mm (see Figs. 2b and 2c), for both the collecting 
and the focusing stages, which are free from spherical aberrations. The latter is an 
essential property, since THz beam spreads very strongly and spherical mirrors are 
thus not suitable. After the collecting mirror, the pump beam is filtered by means of an 
IR- opaque material, whereas the THz beam is focused into the detection crystal -a 
<110> ZnTe sample- together with the optical (IR) probe pulse, thus providing the 
coherent nature of the set-up. Also the focusing lens allows to adjust the waist size of 
the focused THz beam. Finally, the EOS technique (see the following subsection) 
permits to detect the train of THz pulses, by means of the downstream electronic 
instrumentation. 

EOS working-principle. The Electro-Optic Sampling is based on the Pockels effect 
inside electro-optic (EO) crystals, the dual phenomenon of the OR, thus it can be 
described with similar equations. Physically, a linearly polarized optical probe pulse 
E PB (co PB ) -generally obtained by splitting the pump beam- is launched into the 
detection crystal together with the THz wave under the PM condition. However, 
compared to the optical field, THz field is almost constant for several tens of optical 
periods. This means that THz field -which we assume linearly polarized as well- 
provides a change of the refractive index along one axis of the EO crystal, making it 
slightly birefringent. Therefore, while E PB is being propagated through the material, 
its polarization becomes slightly elliptical depending on the THz strength. By 
separating the two optical orthogonal components and sending them to a balanced 
photodetector, THz field time-varying amplitude is recovered by their difference. As 
indicated in Fig. 5a, this process is said "sampling", since the electronics on the 
detection side captures a sample of the instantaneous THz pulse at each laser probe 
pulse so the pulse is simply obtained by bringing together this value sequence. 
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